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We consider the spin-1/2 antiferromagnetic Heisenberg model on three frustrated lattices (the dia- 
mond chain, the dimer-plaquette chain and the two-dimensional square-kagome lattice) with almost 
dispersionless lowest magnon band. Eliminating high-energy degrees of freedom at high magnetic 
fields, we construct low-energy effective Hamiltonians which are much simpler than the initial ones. 
These effective Hamiltonians allow a more extended analytical and numerical analysis. In addi- 
tion to the standard strong-coupling perturbation theory we also use a localized-magnon based 
approach leading to a substantial improvement of the strong-coupling approximation. We perform 
extensive exact diagonalization calculations to check the quality of different effective Hamiltonians 
by comparison with the initial models. Based on the effective-model description we examine the 
low-temperature properties of the considered frustrated quantum Heisenberg antiferromagnets in 
the high-field regime. We also apply our approach to explore thermodynamic properties for a gen- 
eralized diamond spin chain model suitable to describe azurite at high magnetic fields. Interesting 
features of these highly frustrated spin models consist in a steep increase of the entropy at very 
small temperatures and a characteristic extra low-temperature peak in the specific heat. The most 
prominent effect is the existence of a magnetic-field driven Berezinskii-Kosterlitz-Thouless phase 
transition occurring in the two-dimensional model. 

PACS numbers: 75.10.Jm 

Keywords: quantum Heisenberg antifcrromagnet, geometrical frustration, localized magnons, Berczinskii- 

Kostcrlitz-Thoulcss transition, azurite 



I. INTRODUCTION 

The study of frustrated quantum antiferromagnets is 
one of the most active research fields in condensed mat- 
ter physics ii Among them there is a wide class of one-, 
two-, and three-dimensional frustrated quantum Heisen- 
berg antiferromagnets with dispersionless (flat) lowest 
magnon bands. There is currently a great deal of gen- 
eral interest in flat-band systems, since new many-body 
phases can be realized there (see Refs. [2h11| and refer- 
ences therein). 

Interestingly flat-band quantum spin systems^ 3 - admit 
the application of specific methods of classical statisti- 
cal mechanics to study their high-field low-temperature 
behavior. For the application of classical statisti- 
cal mechanics on the quantum systems the concept 
of many-body independent localized-magnon states is 
crucialf 2 - - — >±2r— Typical ground-state features related to 
the localized-magnon states are zero-temperature magne- 
tization plateaus and jumps? 2 - high-field spin-Pcicrls lat- 
tice instabilities^ and a residual ground-state entropy 
at the saturation field4 ' 12 i 13 ' 16 Furthermore, these states 
set an additional low-energy scale that dominates the 
low-temperature thermodynamics in the vicinity of the 
saturation field resulting, e.g., in an extra peak in the 
specific heat at low temperatures.— ^ In two-dimensional 
systems localized-magnon states may lead to a finite- 
temperature order-disorder phase transition of purely 



geometrical origin4^ It is worth mentioning that this 
concept for quantum spin systems is related to Miclkc's 
and Tasaki's flat-band ferromagnctism of the Hubbard 
modeliii^^ 

The previously developed theories for localized- 
magnon spin lattices are valid if the conditions for lo- 
calization of the magnon states are strictly fulfilled (so- 
called ideal geometry). In real- life systems we cannot 
expect that, rather the violation of the localization con- 
dition is typical. Therefore several questions arise: What 
happens when the localization conditions are (slightly) 
violated? Which features of the localized-magnon sce- 
nario survive? Which new effects may appear? 

Although a systematic quantitative theory for such a 
case has not been elaborated until now, it is in order 
to mention here Ref. [l9| considering a distorted frus- 
trated two-leg spin ladder (see also Ref. [20| related to 
a frustrated bilayer) and Refs. [2l| and [22| dealing with a 
distorted diamond spin chain. These studies, however, 
are not based on the localized-magnon picture and use 
a strong-coupling approach, see below. In our recent 
paper— we touch this problem suggesting an heuristic 
ansatz for the partition function of a distorted diamond 
spin chain inspired by localized-magnon calculations. 

The aim of the present paper is to develop a sys- 
tematic treatment of a certain class of localized-magnon 
systems, namely the monomer classji 3 - in the presence 
of small deviations from ideal geometry. We consider 



three different frustrated spin lattices belonging to the 
monomer class, the diamond chain, the dimer-plaquette 
chain (Fig. Q}, and the square-kagome lattice (Fig. [2]). 
We mention that these frustrated quantum antiferromag- 
nets were discussed previously in the literature by various 
authors &M 

Our strategy is to eliminate high-energy degrees of free- 
dom, this way constructing low-energy effective Hamil- 
tonians which are much simpler to treat than the initial 
ones. Among the considered lattices the two-dimensional 
square-kagome lattice is particularly interesting, since 
it shares some properties with the kagomc lattice.— 
Moreover, generally in two dimensions a richer phase 
diagram can be expected. The square-kagome lattice 
also admits a straightforward application of the strong- 
coupling method developed in previous papers, see, e.g., 
Refs. I2lll22ll27l Note, however, that the strong-coupling 
approach is not custom-tailored to the problem at hand, 
since it does not take advantage of the special localized- 
magnon properties. We also want to compare our the- 
oretical findings with available experimental results for 
azurite^ Azurite is known to be the model compound 
for a diamond spin chain (for other compounds with di- 
amond structure, see Ref. [29|). Although its exchange 
parameter set differs from the ideal localized-magnon ge- 
ometry, it is not too far from it . 22 ' 30 

The rest of the paper is organized as follows. First we 
introduce the models and briefly illustrate the localized- 
magnon scenario (Sec. |TTJ) . Then we construct effective 
Hamiltonians (Sec. lIII[) considering separately the strong- 
coupling approach and the localized-magnon based ap- 
proach. In Sec. IHII we compare exact-diagonalization re- 
sults for the full and the corresponding effective mod- 
els to estimate the validity of the effective models. In 
Sec. [TV] we use the constructed low-energy effective mod- 
els to discuss the low-temperature properties of the initial 
frustrated quantum antiferromagnets at high fields. We 
summarize our findings in Sec. |VJ 



II. MODELS. INDEPENDENT LOCALIZED 
MAGNONS 

In the present study we consider the standard spin-1/2 
antiferromagnetic Heisenberg model in a magnetic field 
with the Hamiltonian 
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S Z = J2 S *> J iJ>°- (2- 1 ) 



Here the first sum runs over all nearest-neighbor bonds 
on a lattice, whereas the second one runs over all TV lat- 
tice sites. Note that [S Z ,H] = 0, i.e., the eigenvalues 
of S z are good quantum numbers. The pattern of the 
exchange integrals Jij of the three different frustrated 
lattices considered here is shown in detail in Figs. [JJ and 
[2J On a particular lattice a localized-magnon state 
can be located within a characteristic trapping cell due 
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FIG. 1: (Color online) (a) The diamond chain and (b) the 
dimer-plaquette chain described by Hamiltonian (|2.ip . The 
trapping cells (vertical dimers) for localized magnons are in- 
dicated by bold red lines (J2 bonds). 




FIG 2: (Color online) The square-kagome lattice described by 
Hamiltonian 1)2.11) . The trapping cells (squares) for localized 
magnons are indicated by bold red lines (J2 bonds). 



to destructive quantum interference^ For the diamond 
and dimer-plaquette chains the trapping cells are the 
vertical dimers, see Fig. [IJ for the square-kagome lat- 
tice the trapping cell is a square, see Fig. [3J Owing to 
the localized nature of these states the many-magnon 
states in the subspaces S z = N/2 — 2, . . .TV/2 — n max , 
n max oc TV can be constructed by filling the traps by 
localized magnons. Clearly, all these states are linear 
independent.— Moreover, it can be shown that these 
localized-magnon states have the lowest energy in their 
corresponding ,!5 2 -subspace, if the strength of the antifer- 
romagnetic bonds of the trapping cells Ji exceeds a lower 



bound 2& 

The degeneracy of the localized-magnon states are cal- 
culated via mapping onto spatial configuration of corre- 
sponding hard-core objects on an auxiliary lattice^ ' 12 i 13 
For the frustrated quantum antiferromagnet on the lat- 
tices at hand this hard-core system is a classical gas of 
monomers on a chain or a square latticei 12 ' 13 

It is important to notice that the magnon localization 
occurs due to the specific lattice geometry and hence re- 
quires certain relations between the bonds J^-. For the 
considered traps (single bond or square) this condition is 
fulfilled if an arbitrary bond of the trapping cell and the 
surrounding bonds attached to the two sites of this bond 
form an isosceles triangle, i.e., J\ = J3 = J in Figs.Q]and 
[2j At low temperatures and for magnetic fields h around 
the saturation field h sa .t = hi (hi = J? + J for the dia- 
mond and dimer-plaquette chains and hi = 2 J 2 + J for 
the square-kagome lattice) the contribution of localized 
states is dominating in the partition function^ In the 
present study we deal with the case when the localization 
conditions are violated. To be specific, in what follows we 
consider a violation of the ideal geometry by taking into 
account different values of Ji and J3 but fixing their av- 
erage, i.e., Ji 7^ J3, J\ + J3 = 2 J, sec Figs. [T] and [51 This 
choice is relevant for azurit o 22 ' 30 and it is appropriate to 
illustrate the main point of our considerations. 

For the further treatment of the models we introduce 
a convenient labeling of the lattice sites by a pair of in- 
deces, where the first number enumerates the cells [m = 
1, ... ,Af, M = N/3 for the diamond chain, M = N/4 
for the dimer-plaquette chain, Af = N/6 for the square- 
kagome lattice, N is the number of sites; for the square- 
kagome lattice the cells are enumerated by the vector 
index m = (m x , rn y )] and the second one enumerates the 
position of the site within the cell, see Figs. [1] and [5] 



namely, the fully polarized state \u) and the one-magnon 
state \d). All other sites have fully polarized spins. As 
the magnetic field decreases from very large values, the 
ground state of the cell undergoes a transition between 
the state \u) and the state \d) at the "bare" saturation 
field ho of a cell. The Hamiltonian H is splitted into a 
"main" part H ma x n [the Hamiltonian of all cells and the 
Zecman interaction of all spins with the magnetic field 
ho; for the dimer-plaquette chain it includes also the in- 
teraction along horizontal bonds, J4, see Fig. [TJb)] and 
a perturbation V (the rest of the Hamiltonian H). The 
ground state \ipo) of the Hamiltonian without the con- 
necting bonds, i.e., Ji = for i 7^ 2, at h — ho = 
is 2 -fold degenerate and forms a model space defined 
by the projector P = \<p )(<Po\- For J, 7^ 0, i 7^ 2 and 
h — ho 7^ we construct an effective Hamiltonian H c g 
which acts on the model space only but gives the exact 
ground-state energy. H e g can be found perturbatively 
and is given by^ 3 - - — 



H cS = PHP 



pvy. 
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£0 — £0 



■VP 



(3.1) 



Here \(p a ), a 7^ 0, are excited states of -ffmain- Finally, 
to rewrite the effective Hamiltonian in a more transpar- 
ent form amenable for further analysis it might be con- 
venient to introduce (pseudo)spin operators representing 
the states of each trapping cell. 

Next we present the concrete effective Hamiltonians 
obtained according to the above described procedure for 
the frustrated quantum Hcisenberg antiferromagnets at 
hand. 



Diamond chain 



III. EFFECTIVE HAMILTONIANS 

A. Strong-coupling approach 

The strong-coupling perturbation theory is well estab- 
lished in the theory of quantum spin systems, see, e.g., 
Refs. I2l1l22ll27l Note that it is not necessarily related 
to the existence of localized-magnon states. We begin 
with a brief illustration of the main steps of such an 
analysis. The strong-coupling approach starts from fi- 
nite elements/cells (e.g., dimers or squares) of the spin 
lattice which do not have common sites and have suffi- 
ciently large couplings, here Ji. The number of cells is 
denoted by Af, Af/N < 1. On these cells the spin prob- 
lem is solved analytically. (In the context of localized- 
magnon states these cells will play the role of the trapping 
cells.) The cells are joined via weaker connecting bonds 
Ji, i 7^ 2. The strong-coupling approach is based on the 
assumption that the coupling J2 is the dominant one, i.e., 
Ji/ J2 <C 1, i 7^ 2. Specifically, at high fields considered 
here only a few states of the trapping cell are relevant, 



In the case of the diamond chain the following two 
states of each vertical bond arc taken into account: 
|w) = I tif 2 ) with the energy J 2 /4 - h and \d) = 
(I ti-U) - I I1T2)) /\/2 with the energy -3J 2 /4. Further- 
more, h Q = J 2 - For the projector onto the ground-state 
manifold we have 
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J m — / m ' 



p 
(|t3>(t3|L,'Pm = (|«><«| + |d)<d|) TO .(3.2) 



The set of relevant excited states \if a ), a 7^ 0, is the 
set of the Af2™ states with only one flipped spin | ^3) on 
those sites which connect two neighboring vertical bonds. 
Moreover, e a = Eo + h = £0 + ^2- After introducing 
(pseudo)spin-l/2 operators for each vertical bond 

T z = \ [\u){u\ \d)(d\),T+ = \u)(d\,T- = |d)<u|(3.3) 
Eq. (|3.ip becomes 

Af 

H cS = J2 [J CC^+i + T y m Tl +1 ) hT r ; + C] (3.4) 
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(3.5) 



The effective Hamiltonian in strong-coupling approach 
corresponds to an unfrustrated spin- 1/2 isotropic XY 
chain in a transverse magnetic fields This result coin- 
cides with that one obtained in Ref. [22l 



2. Dimer-plaquette chain 

In the case of the dimer-plaquette chain the rele- 
vant two states of each vertical bond are the same, but 
the projector onto the ground-state manifold obviously 
contains the projector on the up-state for both spins 
on each horizontal bond [J4 bond in Fig. HJb)], i.e., 
P m = V m <g> (| taXta I ® I t4>(t4 |) TO , where V m is de- 
fined in Eq. (|3.2[) . Instead of a flipped spin on the 
site which connects two neighboring vertical bonds be- 
ing an excited state for the diamond chain, we consider 
here two classes of excited states. The first one con- 
tains the singlet state on the horizontal bond, \s) = 
(I T3-W - I -1-314)) /\/2- The energy of this excited state 
is e a — £o + -h — J a- The second one contains the S z = 
component of the triplet state on the horizontal bond, 
I*) = (I T3I4) + I I3T4)) /\/2. The energy of this excited 
state is e a = £0+ J 2 - As a result, we again obtain the one- 
dimensional spin- 1/2 isotropic XY model in a transverse 
field given by Eq. (|3.4p . however, with different parame- 
ters, 
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4 2 16 (J 2 -J 4 ) J 2 ' 
/ii = J 2 + J, J = Jl t ^ ■ (3-6) 



Note that we get J — >• if J 4 — >• as it is expected from 
physical arguments. 



3. Square-kagome lattice 

In the case of the square-kagome lattice we take 
into account the following two states of each square: 



I") = I titatat*) witn the energy J 2 - 2h and \d) = 
(I titata^) - I T1T2I3T4) + I tihisU) - I I1T2T3T4)) /2 
with the energy —J% — h, and ho = 2J 2 . For 
the projector P m in Eq. (|3.2[) we now have P m = 

V m ® (| fsXts I ® I teXte l) m , Pm = (|«X«| + |d)(d|) m) 
see also Fig. [21 Similar to the case of the diamond chain 
the relevant excitations have one flipped spin on those 
sites connecting two neighboring squares either in the 
horizontal or in the vertical direction. Their energy 
is e a = £0 + ho = £0 + 2J 2 . The resulting effective 
Hamiltonian (|3.1[) has the form 

AT 

H CS = JY1 ( T - T « + T ™ T n ) - h E T rn + ^, (3.7) 
(mn) rn=\ 

where the first sum runs over all nearest neighbors on a 
square lattice and 



(■/i - j 3 y 

16J 2 

(Ji - j a f 
8j 2 



h = h — hi — 



2 2 16J 2 



hi = 2J 2 + J,J = 



Ji + Jz 



(3.8) 



Again we get a basic model of quantum magnetism, the 
square-lattice spin-1/2 isotropic XY model in a trans- 
verse magnetic field. 



B. Localized-magnon approach 

As illustrated above, the strong-coupling considera- 
tion provides a low-energy description for the considered 
frustrated quantum Hciscnberg antifcrromagncts at high 
fields provided that the intracell coupling J 2 is much 
larger than all other couplings. This might be a quite 
natural requirement for cell-based (e.g., dimer-based or 
square-based) spin systems but it is not necessary in 
the context of localized-magnon systems. Indeed, the 
localized-magnon picture for the considered spin systems 
emerges when J 4 — J 3 = whereas ( J\ + J 3 ) /2 = J is 
only smaller, but not much smaller than J 2 J^ 

At first glance, we may extend straightforwardly the 
above described scheme introducing another splitting of 
the Hamiltonian H taking advantage of the specific fea- 
tures of the localized-magnon scenario. As the main part 
of the Hamiltonian we take that part of the initial Hamil- 
tonian that corresponds to the ideal geometry, i.e., a 
Hamiltonian with J = ( J\ + J 3 ) /2 instead of J\ and J 3 
at h = h\. The remaining part of the initial Hamiltonian, 
which contains Ji — J and h— hi only, we consider as the 
perturbation. The ground state of the main part is the 
well known set of independent localized-magnon states, 



however, the excited states which are necessary for cal- 
culation of the second term in Eq. (|3.1[) are generally 
unknown. 

We may overcome this difficulty considering as a start- 
ing point instead of H the Hamiltonian 



Eq. dM 



N 



H=THT, ■p = ®% l=1 V, 



(3.9) 



where V m is the projector on the relevant states of the 
trapping cell to. It was defined for the various models 
under consideration in Sec. 1111 Al It is important to note 
that by contrast to the projection operator P used in 
strong-coupling approximation [Eqs. (|3. 1|) and (|3.2|) ]. the 
projection operator V used here does not fix the inter- 
mediate spins connecting the trapping cells thus allowing 
more degrees of freedom. Nevertheless, we have made an 
approximation by reducing the number of states taken 
into account [namely, instead of 4 (16) states of each ver- 
tical bond (square) we consider now only 2 of them, the 
fully polarized state \u) and the localized- magnon state 
\d)}. Moreover, the restriction to states \u) and \d) limits 
the possibility of spreading the cell states over the lat- 
tice. For the reduced set of degrees of freedom it is then 
straightforward to introduce in % again (pscudo)spin op- 
erators. 

The usefulness of this new Hamiltonian T-L is twofold. 
First, it is interesting in its own rights providing an ef- 
fective description of the initial spin model in terms of a 
simpler model with a smaller number of sites, see below. 
Second, a major advantage is that for 7i we can elimi- 
nate the spin variables relating to the sites which connect 
the traps perturbatively with respect to small deviations 
from the ideal geometry (but not with respect to the total 
strength of the connecting bonds) arriving at an effective 
model which is certainly an improvement of the strong- 
coupling one. More specifically, we split % into a main 
part "Hmain (that is the Hamiltonian % with J\ = J% = J 
and h = hi) and a perturbation V = % — H m ain- The 
ground state of H ma in is the same as that of the strong- 
coupling approach |y>o) (although with another value of 
the ground-state energy e , since the connecting bonds 
are present in % m ain) and, hence, also the projector is 
the same P = \ipo)(ipo\. Moreover, for the Hamilto- 
nian 'H m ain ah relevant excited states \tp a )i a 7^ 0, are 
known and therefore an effective Hamiltonian (|3.1[) can 
be worked out. 

Below we present further details for each frustrated 
quantum Heisenberg antiferromagnet under considera- 
tion. 



Diamond chain 



After eliminating irrelevant states of the vertical bond 
in favor of the two relevant ones, |u) and \d), and in- 
troducing pseudospin T-operators (|3.3[) we obtain from 
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That is a spin-1/2 XX Z model with alternating isotropic 
XY bonds in an alternating magnetic z-field on a simple 
chain of 2AT sites, i.e., the unit cell contains two sites. 

To exclude further the spin variables at the sites to, 3, 
to = 1, ...,J\f perturbatively, we consider the main 
Hamiltonian 
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where hi = J2 + J. That is an Ising chain Hamilto- 
nian with known cigenstatcs \(p a ). The set of ground 
states f\ m (|v)| T3))t7»> wnere v is either u or d, has the 
energy e = (— 5 Jif^— J/2)N. Now we consider the per- 
turbation V = H — "Hmain and the set of excited states 
\ip a ) which enter Eq. (|3.1[) . Again we consider the states 
with one flipped spin on the site to, 3. However, now 
(since J / 0) the energy of the excited state depends 
on the states of the neighboring vertical bonds. Namely, 
£ a = £0 + J% — J if for the neighboring sites T z = 1/2, 
£a = £q + J2 if for one of the neighboring sites T z = 1/2 
and for another one T z = —1/2, and e a = Eq + J2 + J if 
for the neighboring sites T z = —1/2. Taking into account 
this circumstance in Eq. (|3.1|) . after straightforward cal- 
culations we arrive at the effective Hamiltonian 



-hT*+C] (3.12) 



with the following parameters 
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We may expand the effective couplings and field with 
respect to J/J2 in Eq. p. 131) : 
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(3.14) 



This result reproduces not only the second-order per- 
turbation theory in 1/ J 2 reported in Refs. [2l| and \2% 
but also the third-order termaSi in I/J2, see Eq. (6.3) of 
Ref . El. 



2. Dimer-plaquette chain 

In the case of the dimer-plaquette chain, the Hamilto- 
nian given in Eq. (|3 .9[) reads: 
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This Hamiltonian corresponds to a spin-1/2 XX Z model 
in a magnetic z-field on a simple chain of 37V sites with 
a unit cell of three sites. 

The starting point for the construction of the TV-site 
effective model is 
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with hi = Ji + J. That is an Ising-Hciscnbcrg chain 
Hamiltonian^ cf. also Eq. (|3. 1 If) . Next, we take 
into account the perturbation V and consider ex- 
cited states \ip a ). In contrast to the strong-coupling 
treatment where the excited states correspond to one 
horizontal bond in the state \s) or |t), now we are 
faced with a more complicated situation, since the 
excited states and their energies depend on the states 
of the neighboring vertical bonds. If both neighboring 
vertical bonds are in the state with T z = 1/2, the 
first excited state | . . . (|«>|*)) m (|«>| t3T4» m+ i • • •) 



has the energy 
excited state 



= £o + J2 — Ji: an d the other 
, .•(l«)l«» TO (l«>Jt3t4» TO+1 ..-> has 
the energy e a = eq + J 2 - Similarly, if both neigh- 
boring vertical bonds have T z = —1/2, the ex- 



cited state | . . . (\d)\s)) m (\d)\ t3T4» m+1 • • •) has 
the energy e a = eo + J 2 — Ji + J, and the ex- 
cited state |...(M)|t)) m (|rf)|t3t4)) m+1 ...) has 
the energy e a — eo + J2 + J- Furthermore, 
if T^ = 1/2 and T z n+l = -1/2, the excited 
state |.-.(| u)Krf)) m (|4lt3t4)) ro +1 ---), \sud) = 
(\s)+x\t))/VT+a?, x = {Ji- ^/Jl + J 2 )/J has the 
energy e a = e + J 2 + J/2 - J 4 /2 - y/ 'J 2 + J 2 /2, and the 
other excited state \ ■ ■ .(\u) \t ud )) m (\d)\ t 3 ti)) m+1 ■ ■ ■), 
= (—x\s) + \t)) j\f\ + x 2 has the energy 



t 



udj 

e a = e + J 2 + J/2 ■■ J4/2 + yJJ 2 + J 2 /2. 
If T^ = -1/2 and T^ +1 = 1/2, the ex- 
cited state |--.(|d) |sd M ))m (l")lt3t4)) m+ i---), 
\s du ) = (|.s) -x\t))/y/l + x 2 has t he energy 
e a = e + J 2 + J/2 - J4/2 - y/Jl + J 2 /2, whereas 
the excited state [ ■ ■ ■ (\d )\tdu)) m (\u)\ t3t4» m+1 ■ • •), 
\t du ) = (x\s) + \t)) j\f\ + x 2 has the energy 
e a = eo + Ji + J/2 - J4/2 + y/J'i + J 2 /2. Taking 
into account all these formulae in Eq. fl3.1[) . after some 
calculations we arrive again at the effective Hamiltonian 
given in Eq. p,12|) . however, now with the following 
parameters: 
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In the limiting case J/J2 —> 0, J / Ji 



, hi = J 2 + J.(3.17) 

-^ this result trans- 



forms into the one obtained using the strong-coupling 
approach, cf. Eq. (|3.6[) . 



3. Square-kagome lattice 

For the squarc-kagomc lattice the Hamiltonian % 
reads 
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.(3.18) 



It corresponds to a spin-1/2 XXZ model on a deco- 
rated square lattice (which is also known as the Lieb 
lattice^), see Fig. [2j The Hamiltonian 'H m ain is given 
by Eq. (|^T5|) with J x = J 3 = J and h = hi = 2 J 2 + J. 
The states with one flipped spin on those sites connect- 
ing two neighboring squares constitute the set of relevant 
excited states. The energy of the excited states depends 
on the states of these two squares. Namely, it acquires 
the value e a = eq + 2 J2 — J if both squares are in the \u) 
state, the value e a — eo + 2 J 2 — J/2 if one of the squares 
is in the \u) state and the other one in the \d) state, 
and the value e a = £0 + 2J 2 if both squares are in the 
\d) state. Taking this into account, we can calculate the 
second term of Eq. (13.11) and arrive at the Hamiltonian 



ing values: 



(Ji - -h) 1 



16J 2 1 



2.1, 



(Ji - -h) / 1 



16J 2 I 1 - JL l-JL- 



2.1 2 



4,/ 2 



h = h — h\ 



(Ji - Jsf 1 




J = Jl ± Js , hi = 2J 2 + J. (3.20) 



Obviously, in the limit J/J 2 — > this result coincides 
with that one obtained within the strong-coupling ap- 
proach, cf. Eq. (pT8j) . 



(mn) 

AT 

-h^T^+TVC. (3.19) 



Here the first sum runs over the neighboring sites of an 
TV-site square lattice and the parameters have the follow- 



C. Exact diagonalization 

In this section we will discuss the region of validity of 
the effective Hamiltonians obtained in the previous sec- 
tions by reducing the dimension of the full Hilbert space 
and by using perturbation expansions. To check the qual- 
ity of the obtained effective Hamiltonians we perform ex- 
tensive exact diagonalization studies for the initial and 
the effective models and compare the results focusing on 
magnetization curves and the temperature dependence of 
the specific heat. 

In Figs. [31 SI and[5l we show the magnetization curve at 
low (including zero) temperatures for the considered dis- 
torted models of TV = 6 (diamond chain, Fig. [3]), TV = 6 
(dimcr-plaquctte chain, Fig. @| as well as TV = 4 and 8 
cells (square-kagome lattice, Fig. [5]). In Figs. [3] and [SJ 



we also show the temperature dependence of the specific 
heat at high fields. Results for the initial full model (|2.1[) 
(thick solid black curves) are compared with correspond- 
ing data for the effective models [(i) strong-coupling 
approach, Hamiltonian H c g, see Eqs. (|3.4[) and (|3.5|l . 
(pO]) and (|3J]l . and (JHHJ) and fl3J§, thin long-dashed 
green curves; (ii) localized-magnon description, Hamilto- 
nian H, see Eqs. ([3~TU]) . ([3"7T5]) , and (|3TS)) . short-dashed 
blue curves; (hi) localized-magnon description with sub- 
sequent perturbation approach, Hamiltonian "H ff, see 
Eqs. { XT!]) and (j3~T3l) . (|XT2j) and ([gTTjl . and (J3~TO|) and 
(|3.20p . dotted red curves]. 

Comparing the full initial model with the effective ones 
we find a quite excellent agreement with the localized- 
magnon based effective models described by H and HeS ■ 
Note that in our plots the deviation from ideal geometry, 
i.e., J\ = J 3 , is up to 20%. We mention that even for 
fairly large deviations, not shown here, the agreement is 
still satisfactorily. Moreover, the perturbative localized- 
magnon description by "H ff almost coincides with non- 
pcrturbative localized-magnon description by % for tem- 
peratures including the full low-temperature maximum 
in c(T). Hence, the former one is the favored approach, 
since it contains much less degrees of freedom and there- 
fore its treatment is simpler. The neglected degrees of 
freedom become relevant only at higher temperatures. 
On the other hand, for J2 = 3 and J2 = 2 considered 
in Figs. [21 SI and[SJ the strong-coupling approach is sig- 
nificantly less accurate. Naturally, it becomes better by 
increasing of J2, see, e.g., the inset in the upper panel of 
Fig. |3l In general, we may conclude that the derived ef- 
fective models approximate the exact results remarkably 
well for a wide range of parameters. 

Some generic features emerging due to the deviation 
from the ideal geometry are already clearly visible from 
the data obtained for the initial full models of small size. 
First we mention that the deviation from ideal geom- 
etry, i.e., J\ ^ J3, does not yield a noticeable change 
of the width of the plateau which precedes the magne- 
tization jump to saturation in the ground state (note 
that the full plateau is not shown in Figs. [3j 01 and 
[51 where data near the upper end of this plateau are 
presented only). However, the deviation from ideal ge- 
ometry has a significant influence on the magnetization 
jump at hi present for J\ = J 3 . Instead of the jump 
there is a (small) finite region, hu < h < hih, where the 
low-temperature magnetization shows steep increase be- 
tween two plateau values. While the strong-coupling ap- 
proach predicts the boundaries of this region hu and h\h 
only qualitatively and underestimates the width of the 
region h\h — hu, the localized-magnon approach gives 
much better results for hu and h\h. Further we no- 
tice that the specific heat c(T) shows a two-peak tem- 
perature profile, typical for localized-magnon systems, 
also for the distorted models. The double-peak struc- 
ture is even present for magnetic fields within the region 
hu < h < h^. All effective Hamiltonians are capable to 
reproduce the low-temperature peak of c(T) related to 
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FIG. 3: (Color online) Comparison of the full model with 
effective models for the distorted diamond chain of Af — 6 
cells: Field dependences of the magnetization (per site) and 
temperature dependences of the specific heat (per site) for 
Ji = 0.85, J2 = 3 (main panels) and J2 = 6 (inset), J3 = 1.15. 
For more explanations, see the main text. 



an extra low-energy scale set by the manifold of almost 
localized-magnon states. Again the localized-magnon ap- 
proach yields much better agreement with the full model 
than the strong-coupling approach. Note that the effec- 
tive model % can reproduce qualitatively even the sec- 
ond high-temperature maximum of the specific heat, see 
Figs. [3] and [5] However, for the description of the low- 
temperature thermodynamics the simpler effective model 
1-L e s is favorable, since both effective localized-magnon 
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FIG. 4: (Color online) Comparison of the full model with 
effective models for the distorted dimer-plaquette chain of 
M = 4 cells: Field dependences of the magnetization (per 
site) for J\ = 0.9, J2 = 2, J3 = 1.1, J4 = 1 (main panel) and 
J 4 = 0.2 (inset) at T = 0.0001. Note that for J 4 = 0.2 (inset) 
all curves practically coincide. For more explanations, see the 
main text. 



based models yield almost perfect agreement with the 
full model. 

Before we will discuss the low-temperature physics of 
the models at hand in more detail in the next section, 
let us summarize the main findings of the present sec- 
tion relevant for the further considerations: We have pre- 
sented three types of effective models which are capable 
to describe frustrated quantum antiferromagnets of the 
monomer universality clasa^ at high fields and low tem- 
peratures. All effective models refer to a reduced Hilbert 
space (only two states for each trapping cell are taken into 
account), the strong-coupling approach assumes J/J2 to 
be small, whereas the localizcd-magnon approach needs a 
more modest requirement: The deviation from the ideal 
geometry should be small. Comparisons of exact diag- 
onalization data illustrate the quality of the suggested 
models. The fairly simple and well investigated spin-1/2 
XXZ model with uniform nearest-neighbor interaction 
in a z-aligned magnetic field on a chain or on a square lat- 
tice, respectively, provides an accurate low-temperature 
description of the frustrated initial spin-1/2 Heisenberg 
models. In the limit of sufficiently strong Ji the even 
simpler spin-1/2 isotropic XY model is adequate. 

We may conclude the above considerations by the 
statement: Instead of a classical hard-core models ad- 
equate to describe the low-energy degrees of freedom for 
ideal geometry we are faced with quantum models in the 
case of non-ideal geometry. However, these quantum 
models are well-known standard ones, they are unfrus- 
tratcd and they allow a straightforward application of 
efficient methods of quantum magnetism such as Bethc 
ansatz, density matrix renormalization group or quan- 
tum Monte Carlo techniques. Moreover, for the effective 
models the particle number is smaller than that of the 
initial models. Finding of effective Hamiltonians is a key 
result of our paper, from which several consequences fol- 
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FIG. 5: (Color online) Comparison of the full model with 
effective models for the distorted square-kagome lattice of 
M — 4 (main panel and left inset) and M — 8 (right in- 
set): Field dependences of the magnetization (per site) and 
temperature dependences of the specific heat (per site) for 
J\ — 0.8, Ji = 2 (main panels and right inset) and Ji — 4 
(left inset), J3 = 1.2. For more explanations, see the main 
text. 



low to be discussed below. 



IV. HIGH-FIELD LOW-TEMPERATURE 
THERMODYNAMICS 

Now we discuss in more detail the properties of the con- 
sidered frustrated quantum antiferromagnets in the high- 
field low-temperature regime using the effective models 
derived in the previous section. For that we use the fact 
that the effective models, the spin-1/2 isotropic XY or 
XXZ Heisenberg model on a chain or on a square lat- 
tice, were extensively studied by different analytical and 
numerical methods in the past. Hence we can use this 
knowledge to discuss the physical properties of the much 
more complicated initial frustrated models in the relevant 
regime. 



A. The one-dimensional case 

The spin-1/2 isotropic XY chain in a transverse field 
is a famous exactly solvable quantum modcl^ There- 
fore, the one-dimensional systems considered here can be 
described analytically even in the thermodynamic limit 
Af — > oo if the strong-coupling approach is adequate, i.e., 
if J2 is sufficiently large. In Fig. [S] we report some ther- 
modynamic quantities for the distorted diamond-chain 
Heisenberg antiferromagnet (|2.1|) with J\ = 0.85, Ji = 6, 
J3 = 1.15 (cf. the inset in the upper panel of Fig- EJ) ob- 
tained by using the effective Hamiltonian H e g given in 
Eqs. d33), (1331) . 

From the upper panel of Fig. [B] it is obvious that 
the ground-state magnetization jump at hi = 7 exist- 
ing for ideal geometry, Ji = J 3 , transforms into the 
smooth magnetization curve of the spin-1/2 isotropic XY 
chain in a transverse magnetic field for the distorted 
model. This smooth magnetization curve has an in- 
finite slope with exponent 1/2 when approaching con- 
tinuously the plateaus at hu = hi = 7 and hih = 
hi + (Ji - J 3 ) 2 /(2J 2 ) = 7.0075. Furthermore, the fi- 
nite residual ground-state entropy per site at the satura- 
tion field, s(T = 0, h = hi) = ln(2)/3 « 0.231, present 
for J 1 = J3, is removed immediately by distortions, i.e., 
s(T = 0, h) = 0. However, by a slight increase of the tem- 
perature the whole manifold of almost localized-magnon 
states becomes accessible for hu < h < hih, thus pro- 
ducing a tremendous entropy enhancement, see the cor- 
responding panel in Fig. [51 The behavior of the specific 
heat, shown in the lower panel of Fig. [6] depends on 
the value of h. For hu < h < hih one has c(T) oc T 
for T — > which corresponds to the gapless spin-liquid 
phase of the effective spin-1/2 isotropic XY chain in a 
transverse magnetic field. Otherwise there is an exponen- 
tial decay of c(T) for T — > (within the gapped phase 
of the effective spin-1/2 isotropic XY chain in a trans- 
verse field) . Correspondingly, the character of maxima in 
c(T) seen in the lower panel of Fig. [6] (representing the 
extra low-temperature maxima of the full model) is dif- 
ferent. It is that one for the spin-1/2 isotropic XY chain 
in the spin- liquid phase when hu < h < hih, but it is 



m(T,h) 




FIG. 6: Magnetization per site m(T,h), entropy per site 
s(T,h), and specific heat per site c(T,h) at high fields and 
low temperatures for the distorted diamond-chain Heisenberg 
antiferromagnet (|2.1|) with J\ = 0.85, J2 = 6, J3 = 1.15 
according to the effective Hamiltonian H c n (|3.4[) . (|3.5[) . 
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the one for (weakly interacting) spins in a field if h is sig- 
nificantly below (above) hu (hi h ). A crossover between 
these types of behavior produces interesting behavior of 
the position and the height of the low-temperature max- 
imum of c(T, h). 

Qualitatively the behavior shown in Fig. [5] remains in 
the case of XX Z (pseudo)spin interactions as well as in 
the two-dimensional case, although quantitative details 
are obviously different [cf., for example, Figs. [TUlandfTTI 
which show results for the (effective) square-lattice XX Z 
model using a slightly different/complementary format of 
representation]. But, clearly, it is much harder to obtain 
analytica l 39 ' 40 or numerical results for XXZ models. 



B. Azurite 

Next, we apply the elaborated effective description to 
discuss some low-temperature properties of natural min- 
eral azurite CU3 (003)2(011)2,— where the magnetiza- 
tion curve is experimentally accessible even beyond the 
saturation fields It is known that the magnetic prop- 
erties of azurite can be explained using the diamond- 
chain Hamiltonian (|2.1j) with the following set of pa- 
rameters: J\ = 15.51K, J 2 = 33K, J 3 = 6.93K with 
h = ou B H, g = 2.06, ^ B « 0.67171K/T, sec Refs. H2 
and [30. Note, however, that also a small exchange inter- 
action J m = 4.62K between the sites m, 3 and m + 1,3 
[see Fig. Q^a)] may be relevant. Neglecting J m the ex- 
change parameters of our effective Hamiltonian T-L s, see 
Eqs. ([332)1 and ([3~13l . are J = 0.845K, J z = 0.287K, and 
the effective magnetic field is h = (1.384H - 44.778)K, 
where H is the experimentally applied magnetic field 
measured in teslas. This effective Hamiltonian is appro- 
priate to describe the low-temperature properties of azu- 
rite above 30T. Note that in Ref. [22| slightly larger values 
of the effective parameters J and J z were obtained (taking 
into account J m = 4.62K) by analyzing the excitations 
above the one-third plateau state and by fitting to the 
values of the upper edge field of the one-third plateau 
and of the saturation field. 

We present results for the azurite parameters for var- 
ious temperatures, including T = 0.08K related to the 
measured magnetization curve,— in Figs.[2|8l and[9l In 
the inset of Fig. [7] we first compare the results of dif- 
ferent approaches for A^ = 6 (i.e., N = 18). Again we 
observe that the results for H, Ti, and Ties agree very 
well, whereas the data for the strong-coupling approach, 
-f/eff, deviate noticeably. Hence we further focus on 
data obtained from our effective Hamiltonian H c h (13.121) . 
(|3~13|) for M = 20, i.e., N = 60, with the parameters 
J = 0.845K, J z = 0.287K, and h = (1.384H - 44.778)K. 
The magnetization curve at T = 0.08K is shown in Fig.[7J 
the field dependences of the specific heat per site c and 
the entropy per site s at T = 0.08K and T = 0.4K 
are shown in Fig. [8j and in Fig. |9] we present the tem- 
perature dependences of the specific heat c and the en- 
tropy s for three values of the magnetic field, H = 31T, 
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FIG. 7: (Color online) Low-temperature (T = 0.08K) mag- 
netization curves for the distorted diamond-chain Heisenberg 
antiferromagnet (|2.1|) with azurite parameters Ji = 15.51K, 
Ji — 33K, Ja = 6.93K (J m = 0), and the gyromagnetic ratio 
g = 2.06. Inset: M = 6, i.e., N = 18; main panel: M = 20, 
i.e., N = 60. 
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FIG. 8: (Color online) Field dependence of the entropy per 
site (thick) and specific heat per site (thin) for the distorted 
diamond-chain Heisenberg antiferromagnet (|2.1[) with azurite 
parameters J\ — 15.51K, J2 = 33K, J3 = 6.93K (J m = 0), 
and the gyromagnetic ratio g — 2.06 at low temperatures, 
T = 0.08K and T = 0.4K. 



32. 5T, and 34T. Remember that the effective model H c s 
covers the low-temperature region, whereas the high- 
temperature maximum in c(T) and the corresponding 
entropy s — > In 2 w 0.693 are not provided by "H ff ■ 

There arc well pronounced features below IK which 
can be attributed to traces of independent localizcd- 
magnon states: A steep increase of the magnetiza- 
tion (Fig. [7|), an enhanced low-temperature entropy 
that reaches ln(2)/3 w 0.231 (Fig. ©, and a low- 
temperature maximum in the specific heat (Fig. [5]). 
The large entropy change caused by variation of the 
magnetic field (see Fig. [5J implies an enhanced mag- 
netocaloric effect as it was noticed already in Refs. [l6| 
and [22J- [Note that recently the magnetocaloric ef- 
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FIG. 9: (Color online) Temperature dependence of the en- 
tropy per site (thick) and specific heat per site (thin) for 
the distorted diamond-chain Heisenberg antiferromagnet (|2.1[) 
with azurite parameters J\ = 15.51K, J2 = 33K, J3 = 6.93K 
(Jm = 0), and the gyromagnetic ratio g — 2.06 at high mag- 
netic fields, H = 31T, 32. 5T, and 34T. 



feet in an XXZ chain has been examined rigorously^ 2 -] 
It is in order to mention here recent experiments on 
the magnetocaloric effect for other (frustrated) quan- 
tum Heisenberg antifcrromagnets4^ Corresponding ex- 
perimental magnetocaloric studies for azurite would be 
of great interest. 



C. The two-dimensional case 

The effective Hamiltonians for the distorted square- 
kagome lattice are well-known spin-1/2 square-lattice 
models: The TV-site square-lattice isotropic XY model 
(|3.7p . the 3TV-site XXZ model on a decorated square 
lattice p,18p . and the TV-site square-lattice XXZ model 
(|3.19p . All models contain a Zeeman term with a mag- 
netic field in z direction. The first and the third models 
have been studied quite extensively in the past (also in 
the context of hard-core bosons; then the magnetization 
corresponds to the particle number and the magnetic field 
to the chemical potential) 42r— 

It is generally known that the classical version of 
the isotropic XY model (without field) undergoes a 
Berezinskii-Kosterlitz-Thouless (BKT) transition^ 2 - at 
T c fa 0.898|J|^ The BKT transition occurs in the quan- 
tum spin-1/2 case too, although, at a lower tempera- 
ture, T c m 0.343|J|i 44 ' 46 ' 48 The quantum model is gaplcss 
with an excitation spectrum that is linear in the momen- 
tum. The specific heat shows T 2 behavior for T — > 0, 
it increases very rapidly around T c , and it exhibits a 
finite peak somewhat above T c . This kind of the low- 
temperature thermodynamics survives for not too large 
z-aligned magnetic field |h| < 2|J| (sec, e.g., Figs. 3 and 
8 in Ref . \5w . Also for the spin-1/2 square-lattice XXZ 
model with dominating isotropic XY interaction in a z- 
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FIG. 10: (Color online) (a) Magnetization per site m(T,h), 
(b) entropy per site s(T,h), and (c) specific heat per 
site c(T, h) at high fields and low temperatures (T = 
0.001, 0.004, 0.01) for the distorted square-kagome Heisenberg 
antiferromagnet (pTT|) with Ji = 0.8, J 2 = 2, J 3 = 1.2. The 
data are obtained from the effective Hamiltonian Ti. c s (|3.19ll . 
(pHD) with TV = 20. 



aligned magnetic field the BKT transition appears.— 

These known results can be translated to the consid- 
ered case, i.e., the distorted square-kagome spin model 
(|2.ip in the high-field low-temperature regime. Again 
we use exact diagonalization for the effective model "H ff 
(|3TT9"|) . (|3~2P|| of TV = 20 sites for the set of parame- 
ters J\ = 0.8, J2 = 2, J3 = 1.2 to analyze the low- 
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FIG. 11: (Color online) (a) Entropy per site s(T,h) and 
(b) specific heat per site c(T, h) at high fields (h — 
4.99,5.01,5.02,5.03) and low temperatures for the distorted 
square-kagome Heisenberg antiferromagnet (|2.ip with J\ — 
0.8, J2 = 2, J3 = 1.2 according to exact diagonalization data 
for the effective Hamiltonian H c s (|37T9]l . ([3~20)) with N - 20. 
Temperature dependence of c(T, h) are shown also for TV = 16 
(thin curves) and M = 10 (very thin curves). 



temperature features of the initial TV-site square-kagome 
model (remember that M — 20 corresponds to TV = 120). 
From the previous section we know that Hog for this set 
of parameters works perfectly well at least up to T = 0.1, 
cf. Fig. O Our results are collected in Figs. [10] and [TT] 
Instead of the former jump of the ground-state magneti- 
zation per site between 1/3 and 1/2 present at h\ = 5 for 
J\ = J3 = 1, the magnetization m(T, h) shows a smooth 
(but steep) increase varying the field h from hu ~ 4.996 
to h\h ~ 5.027, see Fig.fTUTa). Due to the distortion, the 
former 2 -fold degeneracy of the ground states at satu- 
ration field is lifted and, as a result, there is no residual 
entropy at T = 0. Since these 2" energy levels remain 
close to each other, in the field region hu < h < h\h by 
a slight increase of T all these states become accessible 
and the entropy shows clear traces of the former resid- 
ual entropy of size ln(2)/6 « 0.116, see Figs. [TOlb) and 
[TTT a) and cf. also Rcf. [l2|. Concerning the behavior of 
the specific heat of the distorted square-kagome model 
for T — > we use the knowledge for the energy spectrum 
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FIG. 12: (Color online) Sketch of the phase diagram of the 
distorted square-kagome Heisenberg antiferromagnet (Ji = 
0.8, J2 — 2, J3 — 1.2) at high magnetic field (thick dashed 
green line) as it is indicated by the position of the maximum 
T* of the specific heat c(T, h) (thick solid blue curve). Thin 
dashed red curve corresponds to the value of c(T* , h) divided 
by 1000. By filled violet circles the values of hu and h\h are 
indicated. 



of the effective easy-plane XX Z model: The spectrum is 
gaplcss with a linear dispersion of excitations in the field 
region hu < h < h\h and, as a result, the specific heat 
is ex T 2 at low temperatures. Outside that field region 
the spectrum is gapped leading to an exponential decay 
of the specific heat for T — > 0. 

With respect to the BKT transition the low- 
temperature maximum of the specific heat, cf. Figs. llOf c) 
and [TlTb) . deserves a more detailed discussion. For the 
gapped spectrum there is a broader maximum in c(T). 
By contrast, in the gaplcss field region hu < h < h\h, 
the maximum occurs at lower temperatures and it is more 
pronounced, i.e., it becomes peak- like [compare, e.g., the 
curves for h = 4.99 and h = 5.01 in Fig. [TlTb)]. As men- 
tioned above, this well-pronounced maximum is located 
somewhat above the BKT transition point T c . 

Comparing the data for M = 10, 16, and 20, see 
Fig- fTTT b) , one clearly see that the maximum in the c(T) 
curves corresponding to hu < h < hui shows significant 
finite-size effects (in particular the height of the maxi- 
mum increases noticeably with system size TV), whereas 
in the gapped regime the c(T) curves are insensitive to 
the system size also around the maximum. The size de- 
pendence of the height of the maximum can be inter- 
preted as signature of the BKT transition (see also the 
discussion in Ref . |44I ) . 

Using the maximum in the specific heat as an indica- 
tor of the BKT transition in the distorted square-kagome 
Heisenberg antiferromagnet we can construct a sketch of 
the phase diagram of the model in the h - T plane, see 
Fig. [T2] The largest transition temperature T c appears 
for zero effective field h, which corresponds to h m 5.011, 
of the initial model, sec Eq. (|3.20[) . A finite effective field 
|h| > yields a decrease in T c , finally T c becomes zero 
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entering the gapped phase at hu or at h\h- As mentioned 
above, these features appear only at very low tempera- 
tures. By lowering of J2 and/or increasing of \J\ — J$\ 
one can get effective exchange parameters being larger by 
one order of magnitude, see Eq. p.20[) . and correspond- 
ingly enlarged T c . Discussing the field dependence of T c , 
sketched in Fig. [TSJ in terms of the initial model we may 
conclude that in the highly frustrated quantum Hciscn- 
berg antifcrromagnet on the square-kagome lattice with 
deviations from ideal geometry the BKT transition ap- 
pears only for non-zero magnetic fields, i.e., we are faced 
with a magnetic-field driven BKT transition. 



V. CONCLUSIONS 

Motivated by recent experiments on frustrated quan- 
tum Hciscnberg antiferromagnets and recent theories of 
localizcd-magnon systems we have investigated the high- 
field low-temperature regime of various frustrated quan- 
tum spin systems belonging to the monomer class of 
localized-magnon systems. In the case of the ideal ge- 
ometry, i.e., when the localization conditions are strictly 
fulfilled, their low-temperature physics is well described 
by noninteracting (pseudo)spins 1/2 in a field. This de- 
scription is equivalent to the mapping onto hard-core ob- 
jects used in previous papers. Deviations from ideal ge- 
ometry, which are relevant for possible experimental de- 
tection of the features related to localized magnons, set a 
new low-energy scale determined by the distortion of the 
exchange constants. An effective description of the cor- 
responding low-energy physics for the distorted systems 
is given by non- frustrated XX Z models, i.e., collective 
quantum phenomena emerge. New generic features due 
to the deviation form ideal geometry appear for magnetic 
fields in the vicinity of the former saturation field h\ of 
the undistorted systems, namely: (i) the ground-state 
magnetization curve exhibits a smooth (but steep) part 



(instead of perfect jump at hi), (ii) there is a drastic 
enhancement of the entropy at very small but nonzero 
temperatures (instead of nonzero residual ground-state 
entropy at hi), (iii) the specific heat exhibits a power- 
law decay (instead of zero specific heat at h\ and expo- 
nentially vanishing specific heat just in the vicinity of 
hi). It is worth noting that strong variation of entropy 
with varying magnetic field leads to a noticeable magne- 
tocaloric effect at low temperatures around the satura- 
tion field. Moreover, we mention that the characteristic 
extra low-temperature peak in the specific heat survives 
in distorted systems. The most prominent collective phe- 
nomenon is the appearance of a BKT transition driven 
by the magnetic field in the considered two-dimensional 
distorted localizcd-magnon system, the square-kagome 
Heisenbcrg antifcrromagnet. Likely, such a scenario is 
not restricted to this particular model, rather it should be 
also present in other two-dimensional quantum Heisen- 
berg systems built by weekly coupled localized-magnon 
cells. We may conclude that the findings of our paper 
can be useful searching for experimental manifestation 
of localizcd-magnon effects. 
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